Nature of vibrational eigenmodes in topologically disordered solids 
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We use a local projectional analysis method to investigate the eflFect of topological disorder on the 
vibrational dynamics in a model glass simulated by molecular dynamics. Evidence is presented that 
the vibrational eigenmodes in the glass are generically related to the corresponding eigenmodes of 
its crystalline counterpart via disorder-induced level-repelling and hybridization effects. It is argued 
that the effect of topological disorder in the glass on the dynamical matrix can be simulated by 
introducing positional disorder in a crystalline counterpart. 
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I. INTRODUCTION 

Vibrational properties of disordered materials is one 
of the current lively topics of modern condensed-matter 
physicail. Features of the disordered vibrational spec- 
trum, such as the boson peak, the loffe-Regel crossover 
ancLijjJpiaJlipnal localization are being investigated (see 
e.gMiiO). 

Disorder in condensed matter can essentially be clas- 
sified into two basic typesEl: (i) lattice disorder in crys- 
tals, e.g. in substitutional alloys, and (ii) topological 
disorder in, e.g., glasses. Crystalline structures with lat- 
tice disorder can be studied analytically and consider- 
able progress has beeju-achieved in understanding their 
vibrational behaviourou'EI. Structures with topological 
disorder have mainly been in-YflstigatRd .numerically but 
in some aspects analyticallyBml30'E3E3liJtj. Different 
possible ways of describing the origin and nature of vi- 
brational modes, particularly in the low-frequency range, 
in such disordered systems are still under debate. 

The vibrational properties of crystals with lattice dis- 
order can be successfully treated mainly because of the 
existence of a well-defined reference structure, viz. the 
same crystalline structure but without disorder. The 
choice of a reference structure for a topologically dis- 
ordered material is less obvious. In some cases, such 
a choice can be based on the existence of similar local 
order in both the topologically disordered structure and 
its crystalline prototype. For example, a-cristobalite is 
found to.be a good crystalline counterpart for vitreous 
silicaMEZl. The existence of similar structural elements 
in both crystalline and glassy systems is expected to lead 
to similar vibrational dynamics. Similarities in the vibra- 
tional dynamics within these structural units can be re- 
vealed by comparing vibrational eigenmodes of the topo- 
logically disordered structure with those of its crystalline 
counterpart. Such a comparison can elucidate the nature 
of the disordered vibrational eigenmodes and their possi- 
ble generic connection to crystalline modes. Basically, it 



can be expected that vibrational eigenmodes in topologi- 
cally disordered structures have approximately the same 
nature as in crystals with lattice disorderLJEa. Namely, 
they can be regarded as strongly hybridized crystalline 
eigenmodes, shifted in frequency due to disorder-induced 
level-repelling effects. In this paper, we present argu- 
ments, based on a numerical analysis of a representative 
topologically disordered structure, that this conjecture 
holds true, at least, for the structure considered. 



The existence of a crystalline counterpart allows a 
comparison to be made between a glass and its refer- 
ence structure and conclusions to be drawn about the 
effect of topological disorder on the spectrum of vibra- 
tional excitations. Moreover, it is possible to compare 
positional, e.g. quenched thermally-induced, disorder in 
a crystal with topological disorder in the corresponding 
glass in terms of their respective influence on the vibra- 
tional properties. In a computer simulation, positional 
disorder can be introduced into a crystalline structure by 
heating it in the course of a molecular-dynamics run up 
to a certain temperature below the melting point. Then 
the dynamical matrix of an instantaneous configuration 
corresponding to this temperature can be calculated and 
analyzed. Thus it is possible to mimic the main features 
of the vibrational dynamics of topologically disordered 
systems by the vibrational dynamics of their positionally 
disordered crystalline counterparts. In the following, we 
demonstrate that this can be done, at least for some rep- 
resentative structures. 



The rest of the paper is arranged as following. The 
local projectional analysis method is developed and de- 
scribed in Sec. II. The model is described in Sec. |l|. 
Results and conclusions are given in Sec. 1^ and 0, 
spectively. 
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II. LOCAL PROJECTIONAL ANALYSIS 

A comparison of the vibrational dynamics in a glass 
and its counterpart crystalline phase can be based on a 
local projectional analysis. Of course, the equilibrium 
atomic arrangements in the glass and the crystal are es- 
sentially different, and thus we cannot expand the eigen- 
modes in a disordered system in terms of the ccystalline 
eigenmodes, as we can do for disordered latticedZl. Never- 
theless, what we can do in such a situation is to compare 
the local atomic motions in topologically similar con- 
stituent structural units in different frequency ranges in 
the glass and the crystal. 

The local projectional analysis uses the information 
about atomic displacements contained in the vibrational 
eigenmodes {e^^i} with eigenfrequencies {oj^}. Here the 
index i enumerates the atoms and the index ^ labels 
the eigenmodes. Let us introduce the normalized lo- 
cal displacement vector {u^'^} for a structural element 
I: u"pj = e^,j /[pi(lu^)]^^'^ , where j enumerates the atoms 
within the structural unit, and pi{u>^) — I^CjP- 
the case of a crystal, an eigenmode is identified by the 
wavevector k and the dispersion branch For a crys- 
talline structural unit I, a disordered structural unit m, 
a crystalline eigenmode k/3 and a disordered eigenmode 
fi, we can define the squared scalar product: 

A,.„K/3,C.^) = |uguMp , (1) 

where u'^jj, u'l!/' have 3Nu components, being the 
number of atoms within the structural unit. This quan- 
tity depends on the mutual orientation of the crystalline 
and disordered structural units. We are interested in 
finding the crystalline eigenmode which most resembles 
a given disordered mode and, therefore, we choose the 
maximum value (among all possible mutual orientations 
preserving the symmetry) of the squared scalar product, 

^i,m(wk/3, W^) = max{A;,,„(wk/3, Wp)} ■ (2) 

The next step is to average j4;_m(wk/3, w^) over all possible 
disordered and crystalline structural units, and thus cal- 
culate the averaged squared scalar product, A(tiJk;3, w^), 

74(a;k/3,Wp) = ^P/(wk/3)Pm(wM)^',™(^k/3,t^M) ' (3) 
Lrn 

where the local normalization factors, pi{m), take into 
account the weight of the corresponding structural unit 
in the entire mode and P = m?''(^k/3)Pm(w;i). The 
averaged scalar products A{ujkfj, uj^) can be reduced to a 
function of two continuous variables ^(wcryst, Wdis) by an 
interpolation. In this case, ^(wcryst, i^dis) quantifies the 
degree of local similarity between one disordered mode 
with a frequency about wjis and a crystalline mode with a 
frequency about Wd-yst- One way to obtain ^(wcryst) i^dis) 
from A{lo\^p , Lo fj) , which was used in the present work, is 



to compile a histogram for a discrete set of equidistant 
values of c^cryst for an arbitrary set of uj^ = o^dis, i-e. for 
each Wdis, 

A(Wcryst,Wdis) = -^TfT-^ 7 A(wk'/3' , W;,) . (4) 

iV^Wcrystj 

Here k'/3' assume values for which tj^cryst £ < 
^cryst + dhj^ whcrc doj is the width of a histogram bin, 
and iV(ti;cryst) is the number of crystalline states in this 
spectral interval of width du). 

III. STRUCTURAL MODEL 

As a representative example of a topologically disor- 
dered model structure, we consider a single-component 
glass with predominantly icosahedral order (the IC glass) 
constructed by means of molecular-dynamics simiiiation 
with the use of a pair-wise interatomic potentiaO (all 
quantities used in this paper are expressed in Lennard- 
Jones reduced unitsE2l, see Refs. pip^ for more detail). 
The 3A^ X 3A^-dynamical matrix (the number of parti- 
cles is = 16000) for a glassy minimum-energy config- 
uration has been calculated and diagonalized, yielding 
all vibrational eigenvectors, {e^^^}, and eigenfrequencies, 
{wp} for this system. A good crystalline counterpart, ior 
the IC glass is the a phase, a Frank-Kasper crystalc3E3. 

For the local projectional analysis of the type described 
above, we have chosen two Z\A Frank-Kasper polyhe- 
dra (point group Dq^), interpenetrating along the six- 
fold symmetry axes. This structural unit {ZIA2), com- 
prised of 22 atoma; represents a short segment of the -72° 
disclination lind^do. Within the unit cell of the a phase, 
there are 16 interior atoms (centers of ZlA polyhedra) 
of eight partially overlapping Z142 units. In the 16000- 
atom sample of the simulated IC glass, 236 such units 
(some of them partially overlapping) covering 20% of all 
atoms have been identified. These atoms, together with 
their nearest neighbors, represent 52% of all atoms. 

IV. RESULTS 
A. Level-repelling effect 

Fig. 1^ shows the calculated ^(cJcryst, Wdis) for dis- 
ordered eigenmodes from the low-, middle- and high- 
frequency parts of the spectrum (a set of crystalline 
eigenmodes uniformly covering the first Brillouin zone 
with 10'^ points in k-space was used). As seen from this 
figure, ^(wcryst, Wdis), as a function of Wcryst, represents 
a broad distribution centered around Wcryst — i^Syst ; the 
frequency of the crystalline spectrum for which the cys- 
talline modes have the greatest overlap with the modes 
of the disordered system. This suggests that a disor- 
dered mode with a frequency about Wdis is generically 
related to crystalline modes about Lu^yg^; in other words. 
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FIG. 1: Averaged squared scalar product of the atomic dis- 
placements of disordered and crystalline local units renormal- 
ized to the maximum value as a function of the crystalline 
frequency tJcryst for three values of the disordered mode fre- 
quency uidis'- *', 1.12, o; 6.84; x, 16.81. Dashed line: vibra- 
tional density of states (/cryst(tt;cryst) for the a phase. 
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FIG. 2: (a) Reduced vibrational density of states, gdiB{(j->)/uj'^ , 
for the IC glass showing the boson peak (BP), (b) Disorder- 
induced level-repelling effect: the disordered mode frequency 
ojdis versus the characteristic crystalline frequency t^'ryst eval- 
uated using two methods: o, direct; •, background-subtracted 
with LJn — 20; ©, background-subtracted with uiq ~ 10, see 
Ref. Pl. The arrows show the connection between the posi- 
tion of the lowest van Hove singularity (VH) and the boson 
peak (BP) via the local projectional analysis results, (c) Re- 
duced vibrational density of states for the a phase. 



the vibrational motion of the ZIA2 structural units in the 
disordered eigenmode is similar to that of the same struc- 
tural units for the indicated eigenmodes in the a phase. 
Since A{ujcrysti^^dis) as a function of Wcryst is broad and 
asymmetric, particularly at high frequencies, we estimate 
the characteristic frequency of the crystalline modes cor- 
responding to a disordered mode associated with wjis as 



the first moment of this function: 
w^J^gt(wdis) = J ujA{uj,uJdis)dLj/ J A{uj,ujdis)duj . (5) 

In Fig. H(b), ij^cryst' calculated both by the direct method 
and by subtracting a background from ^(wcryst, ^dis)0, 
is plotted versus w^is. 

From Fig. ^ it is evident that the peak-shaped func- 
tion, ^(wcryst, ti^dis), bccomes increasingly broader with 
increasing frequency, Wcryst, meaning that there is an in- 
creasingly weak correlation in character between modes 
in the disordered system and specific similar modes in the 
crystal. Such a large broadening of the peakS|-at high fre- 
quencies has also been found in lattice modelaJ subject to 
large force-constant disorder and, in this sense, it is not 
surprising and reflects the large degree of disorder in the 
IC-glass. A similar effect is known for vibrational plane 
waves (analog of crystalline eigenmodes) propagating in 
glass. The distribution of the weights of different dis- 
ordered eigenmodes contributing to a propagating plane 
wave becomes increasingly broad with increasifttt, plane- 
wave frequency (the strong-scattering regime)Ell3. The 
shape of this distribution can even approach the shape 
of the entire frequency spectrum of a glass. However, at 
low frequencies, and specifically in the boson-peak region, 
Lo-QP — 2.5, the peak width of ^(ajcryst, Wdis) is rather nar- 
row, indicating that there is quite a strong correlation 
between the two types of modes there (this iS|-aiialogous 
to the weak-scattering regime for plane waveaJu). 

The deviation of the curve w^Jyst(ti;dis) from the 

straight-fine bisector, w^ryst = '^dis, reflects the level- 
repelling effect. It is clearly seen that the disordered 
modes from the low-frequency region are pushed down 
in frequency as compared to the related crystalline 
frequency, while the disordered modes from the high- 
frequency part of the spectrum arc shifted upwards with 

respect to Wcryst- 

Why do we describe this phenomenon as a level- 
repelling effect? One can argue alternatively that disor- 
dered sytems are "softer" than their crystalline counter- 
parts, and this results in a smaller sound velocity and in 
excess low-frequency eigenmodes which have nothing to 
do with the level- repelling effect. In the high-frequency 
range, the extra modes could be due to vibrations in local 
structural blocks characterized by relatively large (due to 
statistical fluctuations) spring constants which are also 
not related to the level-repelling effect. These arguments 
are not based on particular asumptions about the topol- 
ogy of the disordered system and therefore should be 
valid for disordered lattices as well. However, we have 
demonstrated the relevance of the level-repelling effects 
in describing vibrational spectra of disordered lattices 
with force-constant disorderQ. This has been done within 
the mean-field approximation (homogeneous disorder), 
ignoring the role of statistical fiuctuations in the force- 
constant distributions. We have explained the appear- 
ance of the low-frequency excess modes (with the VDOS 
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FIG. 3: The reduced vibrational densities of states gi\s(Lo) j iJ^ 
and gcryst(ix')/tJ^, for the IC glass and the a phase, respec- 
tively. The position of the boson peak (BP) in the VDOS of 
the IC glass is marked by the arrow. 



c>c o;^) in terms of the well-known efFectEl of the shift of 
the low-frequency part of the spectrum of quasi-particles 
to lower frequencies due to level-repelling effects. The 
mean-field approach is very accurate in the low-frequency 
regime. This demonstrates the relative unimportance 
of local soft regions (soft configurations) characterized 
by relatively small force constants which are due to sta- 
tistical fiuctuations. These regions, of course, exist but 
the low-frequency eigenmodes are not localized at them 
due to the delocalized properties fAA'^ vibrational states 
spectrum around zero frequencyll3t3. We have also ex- 
amined and confirmed this by means of multifractal anal- 
ysis (to be published elsewhere). This is in contrast to 
the high-frequency part of the VDOS where the mean- 
field approach describes well only the main features and 
tendencies in the spectrum, such as the shift to higher 
frequencies (again due to level-repelling effects) but not 
the band tails containing localized states where statistical 
fluctuations of the force constants are very important. 

Thus, bearing in mind the importance of level-repelling 
effects in determining the spectrum of disordered lat- 
tices, we conjecture that a similar effect can also be im- 
portant in topologically disordered structures. By this 
we mean level-repelling effects relative to the spectra 
of corresponding crystalline counterparts of disordered 
sytems which are caused by positional and topological 
disorder. The crystalline counterparts are characterized 
by the same local order as in the corresponding disor- 
dered structures. In topologically disordered systems, the 
atoms in local structural units are displaced from their 
ideal crystalline positions and the values of the force- 
constants are distributed, giving a situation resembling 
lattice models with force-constant disorder. However, the 
force constants of bonds emanating from a given central 
atom can be correlated with each other (see also Ref. 27). 



This situation is in contrast to lattice models with uncor- 
related force constants. The strength and role of these 
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FIG. 4; Calculated values of the temperature dependence of 
the vibrational contribution to the heat capacity for the IC 
glass (solid line) and the a-phase crystalline approximant (cir- 
cles), plotted as Cv/T'^. The curve for the IC glass calculated 
by omitting all the states with < 4 (i.e. in the boson-peak 
region) is also shown by the dot-dashed curve. The calcu- 
lated Debye values for Cv/T^ for the IC glass (dashed line) 
and cr-phase (dotted line) are also given. 



correlations in the vibrational spectrum are not presently 
known. Another distinctive feature of glassy systems is 
the presence of topological disorder, the role of which in 
vibrational spectra is not yet established. In our current 
analysis, we try to find out to what extent the "lattice 
ideology" (including the level-repelling effect) is applica- 
ble to glassy models. 



B. The boson peak 

Fig. |^(b) can be used to establish the nature of dis- 
ordered modes in different frequency ranges. The low- 
frequency regime is of particular interest for vibrations in 
glasses, because of the so-called boson peak (BP) which 
occurs in the reduced VDOS, g^^Jyio) j lo^ , [see Figs. |^(a) 
and ^ and in the temperature dependence of the (vibra- 
tional) heat capacity Cv, normalized by the Debye 
dependence [the solid line in Fig. |^]. 

The vibrational contribution to the heat capacity per 
particle is given by 



a(T) 



1 / 



T J [exp(fiijj/r) - 1]' 



(6) 



Since all quantities in this paper are expressed in 
Lennard-Jones reduced units, we use the value of /i ~ 
0.02959 referred to argon. The Debye approximation of 
the above equation is given by 



(7) 



where On =rf^D = 2TTh{9N/'i'KVy^^CD is the Debye 
temperaturecl. The Debye velocity of sound cjj is de- 
fined as CD = (1/cf -I- 2/cf)~^^^, where c/ and Ct are the 
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longitudinal and transverse sound velocities respectively, 
which can be estimated from the slopes of the acoustic 
dispersion branchesEil. The Debye velocities for the IC 
glass and the a phase are c^^ ~ 3.43 and c^^^^ ~ 4.74, 
respectively. 

In Fig. ^, the peak in the curve of C^/T^ versus T 
is due to states in the boson-peak range of frequencies 
in the vibrational spectrum. Indeed, if we cut the low- 
frequency part of the spectrum, including the boson peak, 
then the peak in CV/T'^ is significantly suppressed (see 
the dash-dotted line in Fig. |^). The heat capacity of the 
IC-glass can be easily compared with that of the a-phase 
crystal. It can be seen from Fig. ^that the plot of C^/T^ 
for the (T-phase crystal [the circles in Fig. ^ tends to the 
Debye limit (the dotted line) at low temperatures. The 
increase above this level at higher T is due to the influ- 
ence of van Hove singularities in the crystalline VDOS. 
This peak for the cr-phase crystal closely resembles the 
peak for the IC-glass, when the modes in the boson-peak 
region of the VDOS are omitted from the calculations. 
However, the peak in the total curve of C^/T^ for the 
IC glass is significantly shifted to lower frequencies, and 
also approaches a higher value of the Debye limit (the 
dashed line in Fig. Eh because of a lower value of the 
sound velocity. 

Following the arrows in Figs. ||(a)-(c), it appears that 
the vibrational states in the boson-peak region mainly 
correspond to crystalline states in the vicinity of the low- 
est van Hove singularity in the <j phase. Although this 
method of analysis is unavoidably approximate, it does 
seem to confirm the physical picture for the origin of the 
boson peak previously found in the f.c.c. force-constaiit 
disordered latticeQ and in models of disordered silicones. 

Assuming such a scenario for the boson-peak origin, 
we can compare its position, wbp — 2.5, with the posi- 
tion of the lowest van Hove singularity, wvh — 4, and 
thus conclude that ojvh ~ wbp ~ wvh, which is an indi- 
cation of strong disorder. In this model, level repelling of 
low-frequency optic modes to even lower frequencies sig- 
nificantly contributes to the excess mode density in the 
boson-peak region (see below). 



C. Branch- hybridization effects 

The level-repelling effect in glass is accompanied by 
another general effect, namely strong disorder- induced 
hybridization of crystalline eigenmodes. In this effect, 
many crystalline eigenmodes contribute to a particular 
disordered eigenmode (cf. the broad distributions of 
^(wcryst, Wdis) in Fig. gj) . Disordcr-induccd hybridiza- 
tion also involves hybridization between different crys- 
talline vibrational branches. A quantitative characteris- 
tic of such hybridization can be defined via the branch 
hybridization parameter, Wj3{Lj^): 

WfiiuJf,) = ^ A{UJUI3, CJ^)/ ^ ^(c^k/3, ^^m) ■ (8) 
k k/3 
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FIG. 5: Hybridization parameters for acoustic (o) and optic 
(•) branches versus the disordered mode frequency tJdis. The 
boson-peak (BP) frequency is marked by arrow. 



Knowing z«^(a;p), we can say which crystalline branch 
mainly contributes to a particular disordered eigenmode. 
In particular, this is of interest for the low- frequency 
regime, where only the acoustic branches exist in the 
crystal, but where we could also expect an admixture 
of low-lying optic branches in a glass due to hybridiza- 
tion effects. The results of a calculation of the hybridiza- 
tion parameter for acoustic, w^c = S^=i^/3' ^^"^ optic, 
Wopt — X]^=4 ° ''^0' branches are presented in Fig. ||(with 
N-^.c. being the number of atoms in the crystalline unit 
cell; A^u.c. = 30 for the a phase). Here we have averaged 
the hybridization parameter wp over all acoustic and all 
optic branches in order to demonstrate the relative role 
of these types (acoustic or optic) vibrations in disordered 
modes. As seen from Fig. ^ even for the lowest-frequency 
disordered modes, the contribution of optic crystalline 
modes is dominant [wopt(t^dis — 1) — 0.8] which indi- 
cates very strong hybridization effects with the dominant 
number of optic modes. Precisely because of such strong 
hybridization, it is not possible to distinguish between 
the contributions of acoustic branches of different polar- 
izations (transverse or longitudinal). The branches of 
different polarization arc distinct for small wavevectors 
only, but this part of the dispersion is suppressed by the 
dominant contribution from all other wave vectors for 
which the polarization is not well defined. It is also seen 
in Fig. 1^ that the contribution of the acoustic branches 
is maximum at low frequencies and then monotonically 
decays with increasing frequency (whereas the opposite 
behavior is, obviously, found for optic modes), thus in- 
dicating an enhancement of the hybridization with optic 
branches. Such behaviour is qualitatively different from 
the plateau behaviour Jound at low frequencies in the 
f.c.c. disordered latticeQ, where only acoustic modes ex- 
ist, and is due to the existence of the extensive optic 
spectrum in the crystalline counterpart of the IC glass. 
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FIG. 6: Distributions of the dynamical matrix elements, 
scaled by the spectral band width corresponding to the IC 
glass, for the IC glass at T = (solid line), a phase at T = 0.4 
(dashed line), and a phase at T = (bars), (a) Distribution of 
the traces of the diagonal 3x3 blocks. Solid circles correspond 
to matrix elements of Z142 units in the IC glass, (b) Distri- 
bution of the traces (with opposite sign) of the off-diagonal 
3x3 blocks. 



D. Distribution of the matrix elements in the 
dynamical matrix 



Both level-repelling and hybridization effects have 
demonstrated the strong influence of disorder on proper- 
ties of disordered and crystaUine modes, thus indicating 
a strong degree of disorder in the IC glass. A possible 
way to quantify the degree of disorder is to compare the 
distribution of the dynamical-matrix elements in the IC 
glass with that in the a phase (see Fig. The positive- 
definite dynamical matrix consists oi D x D blocks {D 
stands for the dimensionality, D = 3 'm this case) . The el- 
ements of the diagonal blocks are subject to sum-rule cat- 
relations with the elements of the off-diagonal blocksc3. 
Therefore, there are four distinct distributions of the ma- 
trix elements: diagonal and ofUiagonal elements in di- 
agonal and off-diagonal blocksE3. In Figs, ^(a)-(b), we 
have plotted the representative distributions, yOj|^g ^^"^ 

off 
Pdiag 



for the traces of the diagonal elements in diagonal 
and off-diagonal blocks, respectively. The trace of the 
off-diagonal block (taken with opposite sign) for atom i 
interacting with atom j can be associated with the effec- 



tive force constant. 



for interactions 



between these atoms. Similarly the trace of the diagonal 
block for atom i describes the effective force constant, 
i^ii — X^Q^ia^iaj interactions of this atom with all 
other atoms. Two conclusions can be made from a com- 
parison of these distributions for the IC glass and the a 
phase. First, the distributions for the glass can be imag- 
ined as being generically obtained from those for the crys- 
tal by broadening the corresponding (5-functions. Second, 
this broadening is strong enough so that no traces of the 
individual 5-functions remain; this is an indication of a 
strong degree of disorder. 

A natural question is, can we actually reproduce the 
distributions for the glass from the crystalline ones? The 
simplest way to try to do this is to introduce thermal (po- 
sitional) disorder in the a phase by increasing the temper- 
ature, and analyzing the dynamical matrix correspond- 
ing to the instantaneous non-equilibrium state in which 
all the atoms are displaced about their crystalline equi- 
librium positions (positional disorder). We have indeed 
found that, for such a force-constant disordered case, the 
distribution of the dynamical matrix elements (cf. the 
dashed and solid hues in Fig. |6|(a) and (b)) and the VDOS 
for the thermally (positionally)-distorted a phase and the 
IC glasscj are similar. This implies that the main features 
of the vibrational spectrum in a topologically disordered 
system can be essentially reproduced by introducing po- 
sitional disorder in its crystalline counterpart, at least 
for the metallic-type system studied here. The existence 
of topological disorder, as distinguished from positional 
disorder about equilibrium crystalline positions, does not 
seem to play a major role in determining the character of 
the vibrational modes, at least in this case. However, we 
have not directly addressed the situation where there are 
differences in the topological connections of atoms with 
constant force constants. 

We have also found that the distribution of the matrix 
elements for atoms belonging to Z142 units involved in 
the present analysis is similar to the total distribution 
(see Fig. |6|(a)). This confirms the representative charac- 
ter of the Z14:2 units for which we have performed the 
local projectional analysis. 



V. CONCLUSIONS 

To conclude, we have presented evidence that the vi- 
brational modes in a topologically disordered glass are 
generically related to the eigenmodes in the correspond- 
ing crystalline counterpart via disorder-induced level- 
repelling and hybridization effects. In particular, the 
extra states in the low-frequency regime (boson peak) 
appear to correspond to crystalline states in the vicin- 
ity of the lowest van Hove singularities in the crystalline 
VDOS. We have proposed a way of defining the degree 
of disorder (weak or strong) by comparing the distribu- 
tion of the dynamical matrix elements for a topologically 
disordered structure with its crystalline counterpart. It 
appears that the main features of the vibrational spec- 
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trum of a topologically disordered (metallic) glass may 
be reproduced by the vibrational dynamics of its posi- 
tionally (thermally) distorted crystalline counterpart. 
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cryst ' 
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sensitivity of this method can be assessed by changing the 
value of LOQ. Fig. |^(b) shows that the low-ojdis data are the 
least sensitive with respect to the definition of the back- 
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